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Q\ Abstract 

We use a dynamical systems analysis to investigate the future be- 
haviour of Einstein-Aether cosmological models with a scalar field cou- 
pling to the expansion of the aether and a non-interacting perfect fluid. 
13^ The stability of the equilibrium solutions are analysed and the results 

| are compared with the standard inflationary cosmological solutions and 

previously studied cosmological Einstein-Aether models. 

^ 1 Introduction 

Einstein- Aether theory [1, 2 consists of general relativity coupled, at second 
derivative order, to a dynamical timelike unit vector field - the aether. It is 
one of several proposed models of early universe cosmology which incorporate 
I a violation of Lorentz invariance [3] . In this effective field theory approach, the 

t-H aether vector field, u a , and the metric tensor g a b together determine the local 

spacetime structure. 

We shall discuss the late time dynamics of Einstein-Aether cosmological 
models; in particular we explore the impact of Lorentz violation on the infla- 
tionary scenario [1] , which provides one of the simplest ways to describe various 
aspects of the physics of the early universe in standard cosmology. More pre- 
cisely, we study the inflationary scenario in the scalar- vector-tensor theory where 
the vector is constrained to be unit and time like, and investigate whether an 
example of the large class of the inflationary solutions proposed pj are stable 
when spatial curvature perturbations are considered. 
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1.1 Einstein- Aether Cosmology 

In an isotropic and spatially homogeneous Friedmann universe with expansion 
scale factor a(t) and comoving proper time t, the aether field will be aligned 
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with the cosmic frame and is related to the expansion rate of the universe. The 
Einstein equations are generalized by the contribution of an additional stress 
tensor for the aether field. If the universe contains a single self-interacting scalar 
field <p (e.g., a scalar inflation which would dominate in any inflationary epoch), 
with a self interaction potential V that can now be a function of <p and the 
expansion rate 9 = 3a /a = 3H, then the modified stress tensor is [jQ 

Tab = V a 0V fc - (^V c 0V c - V + ev e )g ab + V e (u a u b - g ab ). (1) 

This corresponds to an effective fluid with pressure p and density p of the 
form p = \$ 2 + V - 9V e and p = \^ 2 - V + 9V g + V e , where V{<j>, 9). 
The energy-momentum conservation law or Klein-Gordon eqn. is 

4> + H + = o, (2) 

The augmented Friedmann equation (with the addition of the aether stress 
to the energy density, and where 8irG = 1 = c, and G is the renormalizcd 
gravitational constant) is given by: 

\e 2 = p+U 2 + V-6V e -- 2 , (3) 
3 2 a 2 - 

where k is the curvature parameter, and the Friedmann metric is given by 

ds 2 - dt 2 - a 2 {t){^— + r 2 dd 2 + r 2 sin 2 My 2 }. (4) 
1 — kr z 

[In [2] k was set equal to zero]. The Raychaudhuri eqn. follows from the 
differentiation of the Friedmann eqn. 

1.2 Exponential Potentials 

Exponential potentials of the form Voe _A< ^ arise naturally in various higher di- 
mensional frameworks, such as in Kaluza-Klein theories and supergravity [5], 
and the dynamical properties of the positive exponential potentials leading to 
inflation in the Friedmann-Robertson- Walker (FRW) model have been widely 
studied O [3 |8] . In Einstein- Aether theories there might exist a coupling be- 
tween the scalar field and the aether field through the aether expansion scalar 
9, which requires one to consider more general types of potentials than the 
standard exponential form. In [2] Barrow proposes an ansatz of the form 

n 

V(9,tt>) = y o exp[-A0]+^a r 6l r exp[(r-2)A0/2], (5) 

r=0 

where Vo,A and {a r } are constants, and concludes that there exist power law 
solutions of the Einstein- Aether-generalized combined Friedmann-Klein Gordon 
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system on the form 



0=|lnt, (6) 
a = t B , (7) 
9 = 3^ = 3Bt-\ (8) 

where B is a solution of the polynomial equation 

1 ™ 2 
S+^^ra^SSr 1 ^-. (9) 

r=0 

He presents two potentials for which explicit solutions to ^ can be found, 
namely potentials where ai = Vi and potentials where only 02 is different 
from zero. This corresponds to, respectively, the ordinary exponential potential 
and an exponential potential plus a term quadratic in the expansion. We shall 
consider one of Barrow's examples where the exponential part couples to the 
expansion scalar. 



2 The Model 

We shall study the case: 

V(6,(f>) ^Vne-^ + a iy /Vo8e~^ Xrt > + a 2 6 2 , (10) 

where, for convenience, we have renormalized the constant a\ (the constant ai 
can be absorbed 1\ and will not play an essential role in the dynamical analysis). 
The potential V(6, 4>) can be assumed to be positive definite, but this does not 
imply that the constants ai,a2 necessarily are positive, as can be seen in the 
positive definite potential (e~5 A( ^ — 1(9) 2 where ai is negative and a 2 — ai 2 /4. 

If ai,<22 are small, the potential can be thought of as a perturbation of the 
standard exponential potential. For very large constants 01,02 we can study 
non-perturbative generalizations. 

For the potential ( 10 1 the augmented Friedmann equation ^ becomes 




, m + 1 , o" n* -Tr^TT^ (U) 



where we have normalized the equation with a factor proportional to the square 
of the expansion. The normalized Friedmann equation suggest a suitable set of 
expansion normalized variables: 




3/,. | 

* := ^2(l + 3o 2 ) 0' * ;= V(l + 3o 2 )^^' K ' = (l + 3o 2 )o 2 2 ' (12) 
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assuming that Vo is a positive constant and that a 2 is larger than —1/3. In 
terms of these variables the Friedmann equation assumes the simple form 



1 = * 2 + $ 2 - K, (13) 

and the Raychaudhuri equation (expressed in terms of the deceleration param- 
eter q) becomes 

Using the Raychaudhuri equation one can express the Klein-Gordon equation as 
a hrst order ordinary differential equation completely in terms of the expansion 
normalized variables, and an expansion normalized time: ^ =W~ 1 : 

^- = -(2 - 2* 2 + $ 2 )* + A$ 2 +6(1 - * 2 )$ (15) 

UT 

where the constants a and A are defined through 



_ 3A«i - A /3(l + 3oa) 

a — — ■=, A = Ay . 

2^/2 V 2 

The evolution equation for $ is directly given from the definition of $ and the 
Klein-Gordon and Raychaudhuri equations: 

^ = (1 + 2* 2 - $ 2 - - A*)$. (16) 

UT 



The equations (15 1,(16) constitute an autonomous system of first order dif- 
ferential equations. The invariant set $ = corresponds to a model with a free 
scalar field, and it divides the two-dimensional (2D) state space into a region 
with an expanding universe ($ > 0) and one corresponding to a contracting 
universe (<& < 0). 



The curvature K is determined from the Friedmann equation (13), but an 
auxiliary evolution equation can be derived from (15), (16): dK/dr — [4^ 2 — 



2<j> 2 _ 25^'$]^, showing that K = is an invariant set that also partitions 
the state space into two disjoint regions: a bounded negative curvature region 
(\I/ 2 + <I> 2 < 1), and an unbounded positive curvature region (\I/ 2 + $ 2 > 1). 

Note that a-i does not explicitly appear in the equations; it corresponds to 
a term in the potential proportional to the square of the expansion and can be 
absorbed by a rescaling of the expansion scalar, but the expansion normalized 
system is invariant under such rescalings. 

Defining |A| = y(9 — A 2 + a 2 ), the equilibrium points (pi) of the system are 
given by 
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Equilibria K, q 

pi : * = 0, $ = 0, - 1, 

p 2 ■ * = ±1, $ = 0, 0, 2 

3A+|aA| -Xa 2 +3\aA\ \a\A\ — 9 + 3A 2 — a 2 

77, • \Ir = ! <B = ! — 

P3 ' 9 + a 2 ' (9 + a 2 )a ' ' (9 + a 2 ) 

3X-\aA\ Xa 2 +3\aA\ Xa\A\ + 9 - 3A 2 + a 2 

r>„ • il< = ! S> = ! ! — ! 

Pi ' ' 9 + a 2 ' (9 + a 2 )a ' ' (9 + a 2 ) 

1 ^ -o + v^+S 2A 2 - 6 - a 2 + aV^+S , 

*=A< $ = 2A ' 2A 2 '° 

1 -a-y/WTl -2A 2 + 6 + a 2 + ayfd?Tl 

*=A< $ = 2A ' 2X 2 '° 



The stationary solution p e has $ < for all values of a G R, A e R+, and 
corresponds to a contracting universe. We shall in the following only consider 
expanding solutions with vanishing or negative curvature and will therefore 
ignore this point since it lies outside the region of interest K < 0, $ > 0. Points 
Pi and pi always satisfies these conditions. p\ is always a saddle, p\ saddle or 
source, and p^ always a source. We are most interested in the equilibra ps, 
Pi, and p 5 . The points p 3 , p 4 , and p 5 are contained in this region only for a 
restricted, partially overlapping, range of values in (A, a)-space. 

Equilibrium Point p 3 : 

Range of validity: 

$ P3 > when (A < 3, a > 0), or (A > 3, a < ~\/\ 2 - 9). (17) 
Eigenvalues: 

_ -5|A| 2 + 4A 2 + 3A|aA| ± ^/W+C 
M± ~ 2(9 + a 2 ) ~ 

where 

B = a 4 + 18a 2 + 15a 2 A 2 +81+54A 2 -A 4 a 2 +A 2 a 4 +9A 4 , C = 2A(3A 2 +a 2 +9)|aA|. 

Discussion: The point p 3 is a sink and inflationary when < a and A 2 < 
i(o 2 +6-aVa 2 + 8). 
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Equilibrium Point p^\ 

Range of validity: 

$ P4 > when (a < 0, a 2 > A 2 - 9). (18) 

Eigenvalues: 

_ -5\A\ 2 + 4A 2 + 3\\aA\ ± y/W^C 
yi± ~ 2(9 + a 2 ) ' 

Discussion: The point p\ is a sink when a < and A 2 < ^(a 2 +6— aVa 2 + 8). 
It is inflationary only for the subset of this region where A 2 < \{a 2 + 6 + 
ay/a 2 + 8). 

Equilibrium Point ^5: 

Range of validity: 

$ P5 > always, but K P5 < only when 2A 2 > 6 + a 2 - a\Ja 2 + 8. (19) 
Eigenvalues: 

//± = -1 ± ^\/48 + 22a 2 - 12A 2 + 2a 4 - 2a 2 A 2 - 2a(7 + a 2 - A 2 )\/a 2 + 8. 

Discussion: The point p 5 is a sink when A 2 > | (a 2 + 6 — aVa 2 + 8) . 
Figure [T] shows the relevant properties of the future attractor of the system 
and the stability properties of the equilibrium pints p?,, pi, p§. 
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Sink Properties and Bifurcation Lines Eq. Point P3 
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Figure 1: The first figure shows the properties of the local sink and the curves 
where bifurcations occurs in parameter space. Region I corresponds to a spa- 
tially flat inflationary universe, region II to a flat decelerating universe, and 
region III to a negatively curved universe with a constant expansion rate. The 
next three pictures show the sink (dark)/saddle (medium)/ source (light) prop- 
erties of the individual equilibrium points (only p 3 is ever a source here) . 

3 The Model with Matter: 

If, in addition to the scalar field, there also exists a barotropic perfect fluid with 
linear equation of state, p = (7 — l)p (satisfying | < 7 < 2), we will get extra 
terms in the equations corresponding to the fluid energy density. Assuming 
that there is no transfer of energy between the scalar field and the fluid except 
through gravitation, we get an additional evolution equation for the fluid energy 
density coming from the matter energy conservation equation. A normalized 
energy density is defined through 
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n, *p 

' (l + 3o 2 )0 2 " 

The augmented Friedmami and Raychaudhuri equations become 

1 = + * 2 + $ 2 - K, (20) 

q := -3( ^ + i) = 1(37 - 2)0 + 2* 2 - $ 2 - a*$, (21) 
and the evolution equations become: 

— = -(2 - 2* 2 + $ 2 - 1(37 - 2)^)* + A$ 2 + a(l - * 2 )$, (22) 
ar 2 

^- = (1 + 2* 2 - $ 2 - a*$ - A* + 1(37 - 2)0)$, (23) 

UT 2 

— = (-(37 - 2)(1 - 0) + 4* 2 - 2$ 2 - 2a*$)0. (24) 
dr 

Note that = defines an invariant set of the three-dimensional au- 
tonomous system of first order differential equations. The equilibrium points 
Pi with = from the previous section are also equilibrium points of the ex- 



tended system ( 22 )-( 24 1 . Linearization about the points will result in the same 
eigenvalues as before for the eigendirections contained in the = subspace, 
but will also pick up a new eigenvalue from the influence of the fluid on the 
dynamics. The additional eigenvalues are: 



Pi- /Lt 3 = -(37-2) 

P2(±): /X3 = -3(7 -2) 



P3 ■ fJ-3 



PA- V3= 



3 7 a 2 + 27 7 - 6A 2 - 2\y / a 1 {a 2 + 9 - A 2 ) 
9 + a 2 

3 7 a 2 + 27 7 - 6A 2 + 2A ^a 2 (a 2 + 9 - A 2 ) 



9 + a 2 

P5 ■ = -(37-2) 

Equilibrium points p\ , pi , and ps all gain an eigenvalue that is negative for 
all 7 > §. The new eigenvalue transforms the source(s) to saddle(s) but retains 
p 5 as a sink in the same range as before. How it affects the points p 3 and p^ 
is not immediately obvious, but as we shall see below it does not affect the 
properties of the sinks when 7 is in the range (§,2). 

Equilibrium Point p^: 

The point remains a sink if, in addition to the conditions in the previous section, 
also 

2A 2 + ?A|cL4| < 7 (a 2 + 9). (25) 
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Solving for A 2 we get 

2A 2 < a 2 + 9 7 - |a|Va 2 + 187 - 9 7 2 . 
For 7=| this condition reduces to 

A 2 < ^{a 2 + 6-a\/a 2 + 8), 



which is the same restriction as before. Since the r.h.s. of (25) is increasing 

2 

3 • 



with 7 we have that p% is a sink in the same range as before for all 7 > 2 



Equilibrium Point p^\ 

The third eigenvalue is negative if 



2A 2 - \x\aA\ < 7 (a 2 + 9), (26) 



which implies 



2A 2 < a 2 + 97 + |a| v / a 2 + 187-97 2 . 



For 7=| this inequality also coincides with with one of the previous restrictions 
on A required for p^ to be a sink. Larger 7 gives a weaker restriction and hence 
is pi a sink in the same range as before if we have 7 > | . 

3.1 Equilibrium Points with Vt 7^ 

There are also additional equilibrium points with ^ 0: 
Corresponding to the flat FRW solution is 

P7 ■. # = o, $ = 0, n = 1, 

with eigenvalues 

3 3 

Mi = -7) M2 = 3 7 -2, M3 = -7-3. (27) 

Hence P7 is a saddle in the entire range considered, | < 7 < 2. 

Equilibrium Point pg: 

There is an equilibrium point, p%, representing the Matter Scaling Solution 8 . 
*=i(-«+Va 2 + 97(2- 7 )), * = |[, 

= --ij(97 + « 2 " aV« 2 + 9 7 (2 - 7)) + 1- (28) 
The point p§ nas zer0 curvature and the deceleration parameter is given by 

Qps = -^( 3 7- 2), 
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which is negative for 7 > |. 

Linearization of the system ( 22 ) - ( 24 1 about the equilibrium point pg yields 
the eigenvalues: 



jj(2- 7 )±ivi? + a£, M3 = 3 7 -2, 



(29) 



where 



V = 81 7 2 A 2 - IO87A 2 + 36A 2 - 324 7 3 - 72 7 2 ci 2 + 648 7 2 + 180 7 a 2 - 8a 2 A 2 + 8a 4 
£ = v/a 2 + 9-9(l - 7 )2(36 7 2 - 108 7 + 8A 2 - 8a 2 ). 

The last eigenvalue is negative for all the values of 7 that we consider in this 
paper. The value of //+ can be both positive and negative depending on the 
values of 7, a, and A. But since the terms under the square root in £ always are 
positive and hence £ real, we have that the expression vP + a£ is either positive 
or purely imaginary. The eigenvalue /i_ will therefore always have negative real 
part, and the point p s will always be a saddle. 

4 Notes on More General Potentials 

We also briefly take a look at more general potentials. In particular, we consider 
the potential of the form ([5]) with a single non-zero a r (r an arbitrary integer): 



where a r here is normalized by an appropriate power of Vo- Defining normalized 
variables 5* and $ as before, 




ir-2) 



1 



(30) 




e 



2 



V20 







the normalized Friedmann eqn. becomes 



1 = * 2 + $ 2 + 2(1 - r)a$ 2 - r - K, (31) 

where a = |(3)2a r and we shall define A = ^75^- 

The normalized Raychaudhuri equation then becomes: 




and the evolution cqns. become 

3#' = -(1 + r(r - l)a$ 2 - r )- 1 [2 - 2* 2 + $ 2 + (2 - 3r)(l - r)h&- r 

+ V3\ra (2 - r)M 2 - r )]f + V3A$ 2 - Vi~Xa(r - 2)$ 2 ~'' (32) 

3$' = -(1 + r(r - l)a$ 2 - r )- 1 [V3A* - 1 - 2^ 2 + $ 2 + 2(1 - r)al> 2 " r 

+ V3\ar^$ 2 - r }$. (33) 
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For r = 1 these equations reduce to the equations studied earlier. 

Let us consider equilibrium points with zero curvature. If K = 0, we obtain 



from (31) 



^ = l - & -2(1 -r)&$ 2 " r . (34) 



Setting the right-hand sides of ( 33 1 , ( 32 ) to zero to obtain the equilibrium values 



$, '5 and using (34), assuming $ ^ 0, we find after some algebra that $ satisfies 



the polynomial eqn. 



\2 

<l 2 = 1 - 2(1 - r)h&- r (1 + r5$ 2 - r ) 2 , (35) 

3 



whence we obtain 



* = ^-{l + rm 2 - r ). (36) 
V3 



It can easily be checked that these values for iff do indeed satisfy the zero- 



curvature condition ( 34 ) and are equilibrium values of the system ( 33 ) . For van 



ishing a, we obtain the usual zero-curvature inflationary power-law equilibrium 



solution P in standard exponential potential cosmology (4> = y 1 — ^A 2 , $ = 
-4|), which is asymptotically stable to the future. Linearization of the system 



<|32j) , |33j) is more complicated with a general power r, but an approximate 
linearization for small a around the zero-curvature solution is possibl^] 

5 Discussion 

We have analyzed the dynamical evolution and stability of inflationary solutions 
of homogeneous and isotropic Einstein- Aether cosmologies containing a scalar 
field, under the assumption that the scalar field interacts with itself and the 
aether through a potential given by ( 10 ). The potential is of the general form ([5| 



proposed in [5] and more general than the examples explicitly studied therein. 

We find that the scalar field-aether interaction term only slightly affects the 
stable solution of the system for small values of the normalized scalar field self 
interaction coupling A, even when the normalized scalar field-eather coupling 
a is large. When large scalar field self interactions are considered there is a 
qualitative change in the future stable equilibrium of the system compared to 
the ordinary exponential potential case. For negative values of the coupling a 
there exists a region (region II, Fig. [I]) where the stable equilibrium is spatially 
flat but non-inflationary, a situation that in the normal case requires fine tuning 
of A. Even larger scalar field self interactions destabilizes the flat equilibrium 
and drives the state towards a spatially curved equilibrium, like in the normal 
exponential potential case. 

When we also include a matter source term in the form of a perfect fluid 
with linear equation of state that does not couple directly to the other fields, 



A work in progress by M. Stevens. 
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we find equilibrium states corresponding to the usual FRW model and matter 
scaling solution. The precise value of the normalized field variables of the latter 
equilibrium depend on the coupling a, but the qualitative properties like its 
curvature and deceleration parameter do not. The equilibrium point is a saddle 
when curvature perturbations are considered, like in the normal case, but even 
within the flat models it can remain a saddle for some values of A and a, unlike 
the case with an ordinary exponential potential where the matter scaling solu- 
tion is a late time attractor. The matter source term do not qualitatively alter 
the stability of the sinks found earlier, they all obtain a stable manifold of one 
dimension larger than previously and thus remain sinks. 

Also potentials with a coupling between the exponential part and the expan- 
sion of general order are susceptible to a dynamical systems analysis in the scale 
invariant variables. The problem of finding a scale invariant solution with zero 
curvature can be reduced to solving a polynomial equation, and the solutions 
all reduce to the standard inflationary power-law solution when the coupling 
becomes arbitrarily small. 
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